An inclusive search for the Higgs boson in the four-lepton final state at CDF 
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An inclusive search for the standard model Higgs boson using the four-lepton final state in proton- 
antiproton collisions produced by the Tevatron at ^/s — 1.96 TeV is conducted. The data are 
recorded by the CDF If detector and correspond to an integrated luminosity of 9.7 fb"'^. Three 
distinct Higgs decay modes, namely ZZ, WW, and rr, are simultaneously probed. Nine potential 
signal events are selected and found to be consistent with the background expectation. We set 
a 95% credibility limit on the production cross section times the branching ratio and subsequent 
decay to the four lepton final state for hypothetical Higgs boson masses between 120GeV/c^ and 
300 GeV/c^ 



The vector gauge bosons mediating the weak force, the 
W and Z, are massive. Within the standard model of 
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particle physics (SM) their masses arise through sponta- 
neous symmetry breaking (SSB) of the electroweak sym- 
metry [1-3] through the introduction of a scalar field that 
maintains the gauge invariance of the theory 043] ■ Fun- 
damental fermions may acquire mass through Yukawa 
couplings with this field. Quantization of the field ne- 
cessitates the existence of an associated spin-0 particle, 
known as the Higgs boson, the discovery of which would 
confirm the SM mechanism for electroweak SSB. 

Direct searches at LEP combined with recent 
search results from the Tevatron and LHC experi- 
ments [13, U], exclude all potential SM Higgs masses 
outside the ranges 116.6 - 119.4 GeV/c^ and 122.1 - 
127GeV/c^. Although signal processes leading to four- 
lepton final states are not the most dominant within these 
mass ranges the inclusion of additional channels improves 
overall search sensitivity. Moreover, using the new tech- 
niques presented here, the search is simultaneously sensi- 
tive to multiple production and decay modes. Also based 
on this inclusive approach, information regarding Higgs 
boson couplings to both fermions and bosons can be ex- 
tracted, which could be useful to probe the SM or beyond 
SM nature of a hypothetical signal. 

Searches for the SM Higgs boson to four leptons 
through the decay ZZ^*^ have been published by AT- 
LAS [l^ and CMS [3] using proton-proton collision data 
at a center of mass energy of 7 TeV. In this Letter a 
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search for the Higgs boson produced inchisively and de- 
caying into final states containing either four electrons 
(4e), four muons (4/i) or two electrons and two muons 
(2e2/i) is reported. Data from proton-antiproton colli- 
sions at center of mass energy 1.96 TeV, collected with 
the CDF II detector IJ] at the Tevatron and correspond- 
ing to 9.7 fb~^ of integrated luminosity, are used. Higgs 
boson decays to Z boson pairs [H — >■ ZZ) are the domi- 
nant contribution for most of the considered Higgs boson 
mass range. Owing to the full reconstruction of the fi- 
nal state, the ZZ channel provides good search sensitivity 
in spite of the small potential signal yields because the 
resonance structure of the signal can be exploited to dis- 
tinguish it from non-resonant backgrounds. Along with 
the direct Higgs production mechanisms of gluon-fusion 
{ggH) and fusion of vector bosons emitted from the in- 
coming partons (VBF), the search is sensitive to asso- 
ciated Higgs boson production processes (VH). In the 
case of ZH production, additional potential four-lepton 
event contributions originate from Higgs boson decays to 
W boson pairs {ZH — > ZWW — > IHvlv) and r-lepton 
pairs {ZH — )• Ztt — > lllvvlvv). The detection of four 
leptons provides one of the cleanest signatures available 
at a hadron collider. Because of the small probability 
associated with mimicking the signature of isolated lep- 
ton candidates in the detector, background contributions 
from ubiquitous multi-jet production processes are neg- 
ligible in this final state. 

Separation of potential signal and background event 
contributions is obtained from the observed four-lepton 
invariant mass {mu) and missing transverse energy spec- 
tra [sif . Sensitivity is significantly enhanced through 
the inclusion of information since H — > WW and 
H ^ TT decays to leptons result in final state neutrinos, 
while none are expected in background events. The inclu- 
sion of I^j, significantly improves the search sensitivity 
at Higgs masses below 150 GeV/c^, and represents the 
main advancement in analysis technique over previously 
published searches in this final state. 

Our search is performed using the CDF II detector con- 
sisting of a solenoidal spectrometer with a silicon tracker 
and an open-cell drift chamber (COT) surrounded by 
calorimeters and muon detectors [l3 |. The geometry is 
characterized using the azimuthal angle 4> and the pseu- 
dorapidity rj. Transverse energy, Et, is defined as E sin 9, 
where E is the energy of an electromagnetic or hadronic 
calorimeter energy cluster. Transverse momentum, px, is 
the track momentum component transverse to the beam 
line. 

Electron candidates are identified by matching a cen- 
tral or forward track to energy deposited within the 
calorimeter. Muon candidates are formed from charged 
particle tracks matched to minimum ionizing energy de- 
position in the calorimeter, which may or may not be 
matched to track segments {stubs) in the muon chambers 
situated behind the calorimeters. Lepton reconstruction 



algorithms are well validated and described in detail else- 
where 1J| . Taus are included in this search only if they 
decay to electrons or muons. 

Candidate leptons are separated into various cate- 
gories: electrons, muons, and isolated tracks that project 
to detector regions with insufficient calorimeter cover- 
age for energy measurements. Electron candidates are 
distinguished by whether they are found in the central 
or forward calorimeters (|ry| > 1.1) where only silicon 
tracking information is available. The electron selec- 
tion relies on track quality, track-calorimeter matching, 
calorimeter energy, calorimeter profile shape, and isola- 
tion information. Most muon candidates rely on direct 
detection in the muon chambers, which are distinguished 
by their acceptance in pseudorapidity: central muon de- 
tectors (1 77 1 < 0.6), central muon extension detectors 
(0.6 < I?/! < 1.0), and the intermediate muon detector 
(1.0 < I77I < 1.5). Remaining muon candidates rely on 
track matches to energy deposits consistent with a mini- 
mum ionizing charged particle in the central and forward 
electromagnetic calorimeters respectively, failing to have 
an associated stub in the muon sub-detectors. Isolated 
tracks pointing to uninstrumented regions of the detec- 
tor and satisfying high-quality requirements are recon- 
structed as leptons of unidentified flavor. All leptons are 
required to be isolated by imposing the condition that 
the sum of the transverse energy of the calorimeter tow- 
ers m a cone of AR = v/(A0)2 + (Ary)2 = 0.4 around the 
lepton is less than 10% of the electron Et (muon px). A 
similar isolation requirement is applied based on the re- 
constructed tracks in a AR — 0.4 cone around the lepton 
candidate track. 

The probability that a jet will be misidentified as a 
lepton is measured using jet-enriched samples and is cor- 
rected for the contributions of leptons from W and Z 
boson decays. An average value is obtained from sam- 
ples collected with different single jet Et threshold re- 
quirements (20, 50, 70, and 100 GeV) and an uncer- 
tainty is assigned based on the spread within the indi- 
vidual measurements. The range of measured misiden- 
tification probabilities for the lepton categories, which 
vary according to Et or pT, is 0.5% — 3% (central elec- 
trons), 2% - 6% (forward electrons), 0.5% - 4% (central 
muons), 0.5% — 2% (extension muons), 0.5% — 2% (in- 
termediate muons), 0.5% — 6% (calorimeter only muons), 
and 0.5% — 3% (isolated tracks). 

Events for the analysis are collected using online event 
selection (trigger) requirements corresponding to the 
presence of a single high-iS^ electron or high-p-r muon. 
The electron trigger requires an electromagnetic energy 
cluster in the central calorimeter with Et > 18 GeV ge- 
ometrically matching the direction of a charged particle 
reconstructed in the COT with pT > 8 GeV/c. Muon 
triggers are based on track segments in the muon cham- 
bers geometrically matching the direction of a charged 
particle reconstructed in the COT with pT > 18 GeV/c. 
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Trigger efHciencies are measured using samples of lep- 
tonic Z decays 2M ■ ensure a uniform trigger efficiency 
over the lepton momentum spectra, an offline selection is 
applied to the lepton matched to the trigger object, re- 
quiring Et > 20 GeV (pt > 20 GeV/c) for electrons 
(muons). 

Additional charged electrons (muons) are required to 
have Et > 10 GeV (pr > 10 GeV/c). Exactly four 
leptons are required, each separated from the others by 
AR > 0.1. This analysis evolved from a CDF measure- 
ment of the ZZ production cross section in the four lepton 
final state [l^, where constraints on the invariant mass 
of opposite-sign same-flavor dilepton pairs were imposed 
in order to explicitly reconstruct Z bosons. For Higgs bo- 
son masses smaller than 180GeV/c^, at least one of the 
Z bosons is off-shell and such requirements on the mass 
become inefflcient. We therefore require dilepton pair 
masses to be between 20 and 140 GeV/c^. In the final 
state where all leptons are of the same flavor, opposite- 
sign pairings are assigned based on the the smallest de- 
viation of the reconstructed masses from the known Z 
boson mass. Because the backgrounds are modest, loose 
constraints on the mass improve the sensitivity of the 
search. Higgs boson production can also involve hadronic 
jets originating from the associated vector boson in WH 
or ZH production, the additional quarks in the forward 
detector from VBF production, or initial state gluon ra- 
diation. Therefore, we place no restriction on the number 
of reconstructed jets contained within each event. 

The selected events consist primarily of the back- 
ground from non-resonant diboson production of Z* /Z- 
boson pairs (ZZ). Smaller contributions originate from 
non-fully leptonic ZZ decays and Z^ production, which 
lead to events with two or three real reconstructed lep- 
tons and one or two falsely reconstructed leptons associ- 
ated with the photon or additional hadronic jet activity. 
The background from top-quark pair production is found 
to be negligible (<0.01 events). 

The acceptances, efficiencies and kinematic proper- 
ties of the signal and background processes are deter- 
mined from simulation, pythia [31 is used to model all 
Higgs boson production processes and the non-resonant 
ZZ background. The modeling of Zj is based on the sim- 
ulation framework developed by Bauer and Berger [l7| . 
For all samples CTEQ5L parton distribution functions 
(PDFs) are used to model the momentum distribution of 
the initial-state partons [l3|. The cross sections for each 
process are normalized to next-to-next-to-leading order 
(NNLO) calculations with logarithmic resummation for 
ggH [Illi^, NNLO for VH [2ll-ii]^nd next-to-leading 
orderTNLO) calculations for VBF [2l|, [11] , [111 , and 
Zj 

The response of the CDF II detector is modeled with 
a GEANT-based simulation [27| . Efficiency corrections for 
the simulated CDF II detector response for leptons and 
photon conversions were determined using independent 



data samples. 

The overall normalization of the estimated background 
from false leptons is derived using a data-driven ap- 
proach while the shape of the distribution is derived from 
MC simulation. The total contribution from false lepton 
events is estimated using events with two or three leptons 
and additional jets that may be falsely reconstructed as 
leptons weighted by the measured jet-to-lepton misiden- 
tiflcation probabilities. Since the number of available 
three-lepton events in data is small, we model kinematic 
distributions using a weighted sum of those obtained from 
the simulated ZZ and Zj MC samples. 

The kinematic distributions of ZZ events with falsely 
reconstructed leptons is assumed to be the same as that 
for correctly reconstructed ZZ events. We apply the same 
procedure to the simulated Z"f sample to obtain the in- 
variant mass distribution for events from this process, 
which are found to be well modeled by a Landau func- 
tion. We model in Z7 events using the distribution 
from simulated ZZ events with falsely reconstructed lep- 
tons, which was found to agree. 

In the ranges 50 < m^i < 600GeV/c^ and < 
l^rp < 200 GeV, we estimate background contributions of 
10.59± 1.34 ZZ events and 0.39±0.19 events with falsely 
reconstructed leptons. For a SM Higgs boson with a 
mass of 125 GeV/c^, contributions of 0.053 (ggH), 0.003 
(VBF), 0.006 (WH), and 0.089 (ZH) events, yielding a 
total of 0.15 ± 0.01 events, are expected. The indicated 
uncertainties include statistical and systematic contribu- 
tions added in quadrature. Systematic uncertainties are 
described below. 

We observe a total of nine events, which is consistent 
with the rate expected from background sources only. 
The four-lepton invariant mass and I^j- distributions are 
shown in Fig. [lla) and Fig. [ijb) respectively, both over- 
laid with expected contributions from backgrounds and 
Higgs boson production for a mass of 125 GeV/c^, shown 
separately for each production process. 

A variety of possible systematic effects were consid- 
ered including those that affect the normalization and 
the shape of the kinematic distributions. The dominant 
systematic uncertainties are those on the theory predic- 
tions for the cross sections of signal and background pro- 
cesses. Systematic uncertainties associated with the MC 
affect Higgs signal, ZZ, and Z7 acceptances determined 
from the simulated event samples. 

Uncertainties originating from lepton selection and 
trigger efficiency measurements are propagated through 
the acceptance calculation, leading to uncertainties of 
3.6% and 0.5%, on the predicted signal and background 
event yields, respectively. In addition, all signal and 
background estimates obtained from simulation have an 
additional 5.9% uncertainty originating from the mea- 
surement of the luminosity [2^ . 

The gg ^ H cross-section has been computed at 
NNLO and next-to-next-to-leading log (NNLL) preci- 
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sion varying the renormalization and factorization scales 
and implementing different correlated combinations of 
the MSTW2008 error PDFs together with allowed val- 
ues of the strong couphng constant (a^) 2^, S^, yielding 
systematic uncertainties of 7.0% and 7.7%, respectively. 
Uncertainties on VBF and associated Higgs boson pro- 
duction, which account for about a quarter of the total 
Higgs boson events in our sample, are of 5% and 10%, 
respectively (sij. A 3% uncertainty was assigned on the 
branching fraction for H — > ZZ and H WW, which 
are 100% correlated, as well as a 3% uncorrelated un- 
certainty on H ^ TT [3l|. The PYTHIA MC simulation 
for ZZ production that is used for determinin g a ccep- 
tances is based on an expansion at LO; mcfm [321 '^^^ 
used to estimate the potential difference in the accep- 
tance arising from a full NLO simulation, estimated to 
be ±2.5% and assigned as a systematic uncertainty. We 
assign a 10% uncertainty on the ZZ cross-section based 
on the difference between LO and NLO fssi predictions. 
The Zj component was parametrized using a Landau 
distribution, where an uncorrelated 50% uncertainty was 
assigned to the yield prediction in each bin to fully cover 
the possible effects of mis-modeling. Misindentification 
probabilities were measured in several jet samples and 
the maximum spread between these measurements was 
assigned as a systematic uncertainty on the background 
estimation. Propagated through to the acceptance, this 
technique results in a 50% variation in the estimate of 
the background from fake leptons. The missing trans- 
verse energy is scaled up and down by 20% to account 
for potential MC mis-modeling, and the modified shapes 
are taken as shape uncertainties on the MC I^j^ distribu- 
tions used as inputs to the limit calculation. 

Only three bins (0-15, 15-45, and 45-200 GeV) 
are used in the limit setting procedure. The varied spac- 
ing was chosen as a compromise between search sensitiv- 
ity and ensuring sufficient number of simulated events in 
each bin. 

As a cross-check the distribution of the number of jets 
per event in data was compared to the expectation from 
the ZZ MC simulation and found to be consistent. 

Upper limits at the 95% credibility level (C.L.) are 
set on the Higgs boson production cross section, aH, as a 
function of to// . A Bayesian technique [s^l was employed, 
where the posterior probability density was constructed 
from the joint Poisson probability of observing the data 
in each bin of the mn—Ifrp space, integrating over the un- 
certainties of the normalization parameters using Gaus- 
sian priors. A non-negative constant prior in the signal 
rate was assumed. The expected limit and associated 
one and two sigma bands are shown along with the ob- 
served limit in TableUand Fig.[2j For a Higgs boson with 
niH = 125 GeV/c^ we expect a sensitivity of 26.5 x asM 
at 95% C.L. while we observe a limit of 29.3 x asM at 
95% C.L. . The results obtained in the investigated mass 
range are consistent with the observed event excess near 




M„ [GeV/cl 




20 40 60 80 100 120 140 160 180 200 
(b) ^T[GeV] 



FIG. 1: Distributions of the four-lepton invariant mass (a) 
and the missing transverse energy (b) in data. The combined 
estimated contribution from non-resonant ZZ production and 
fake leptons, denoted as background, is overlaid. The poten- 
tial contributions of the different Higgs production processes 
for niH ~ 125 GeV/c'^ are stacked and also overlaid. 



325 GeV/c^ in the CDF search for high-mass resonances 
decaying to ZZ [35]. The analysis reported here is per- 
formed using standard CDF tracking algorithms while 
[ssf uses an alternative reconstruction. 

In 9.7 fb^^ of integrated luminosity collected at the 
Tevatron, no evidence for a Higgs boson signal is found 
in the mass range explored. Upper limits on the Higgs bo- 
son production cross section in the inclusive four-lepton 
final state relative to the SM expectation are obtained. 
Maximal search sensitivity is obtained not only in the 
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TABLE I: Expected and observed upper limits on the Higgs boson production relative to the SM prediction for Higgs particle 
masses from 120 GeV/c^ to 300 GeV/c^ 



H-»4£ 120 13n 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300 

Expoctod/CTq;^,f 38.0 18.3 11.7 9.4 16.0 25.1 18.5 9.8 10.6 12.0 15.7 16.6 18.9 20.5 21.1 23.2 23.5 28.0 30.5 
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FIG. 2: Expected and observed upper limits on the Higgs 
boson production relative to the SM prediction as a function 
of Higgs boson mass. 

high mass region where both Z bosons are produced on- 
shell but also in the lower mass region around 150 GeV/c^ 
where additional signal contributions from ZH — >■ ZWW 
and ZH Ztt improve the sensitivity by 15%. 
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